Experimental investigation of planar ion traps 
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Chiaverini et al. [Quant. Inf. Comput. 5, 419 (2005)] recently suggested a linear Paul trap geometry for ion 
trap quantum computation that places all of the electrodes in a plane. Such planar ion traps are compatible with 
modern semiconductor fabrication techniques and can be scaled to make compact, many zone traps. In this 
paper we present an experimental realization of planar ion traps using electrodes on a printed circuit board to 
trap linear chains of tens of 0.44 jUm diameter charged particles in a vacuum of 15 Pa (10"^ torr). With these 
traps we address concerns about the low trap depth of planar ion traps and develop control electrode layouts 
for moving ions between trap zones without facing some of the technical difficulties involved in an atomic ion 
trap experiment. Specifically, we use a trap with 36 zones (77 electrodes) arranged in a cross to demonstrate 
loading from a traditional four rod linear Paul trap, linear ion movement, splitting and joining of ion chains, 
and movement of ions through intersections. We further propose an additional DC biased electrode above the 
trap which increases the trap depth dramatically, and a novel planar ion trap geometry that generates a two 
dimensional lattice of point Paul traps. 



1. INTRODUCTION 

In recent years, the quantum computing community has 
demonstrated the basic building blocks of a scalable ion trap 
quantum computer 1 1, 2, 3, 4]. The design uses the electronic 
states of ions trapped in a radiofrequency (RF) Paul trap as 
qubits and accomplishes logic gates and state readout by laser- 
ion interactions. Ions are shuttled between zones in an array 
of ion traps to perform two qubit gates between arbitrary pairs 
of qubits. Recent experiments have accomplished state read- 
out 1 5], one and two qubit gates 00], and ion shuttling in 
straight lines and through tees OTBl- Scaling this architec- 
ture up to the many thousands of qubits necessary for a use- 
ful computation, however, will involve significant physics and 
engineering challenges. 

One of these challenges is to build a many zone ion trap 
array capable of holding several thousands of ions in memory 
zones and moving arbitrary pairs of them together in interac- 
tion zones. The ion traps used for quantum computation are 
based on the linear RF Paul trap [10^, II, I2j 1^]- Current 
experiments typically use gold electrodes deposited on two or 
more alumina substrates with a geometry similar to that shown 
in Fig. [Ha) 1 14, 15]. Ions are confined radially to the null axis 
of an RF quadrupole electric field, while axial confinement 
and ion movement operations are accomplished by DC elec- 
tric fields. While these multi-level traps can be adapted to 
microfabrication techniques 1 16], it is not clear whether they 
can be scaled to many zone traps because they require slots 
through the trap structure. This will make trap topologies that 
include loops difficult, as there will be islands of electrodes 
that will have to be mechanically supported and electrically 
connected. 

Chiaverini et al. 1 17] proposed using a planar RF Paul trap 
geometry for ion trap quantum computing which is easy to 
scale up to many zone traps and amenable to modern micro- 
fabrication techniques. The electrodes all lie in a plane and 
ions are trapped above the plane of the electrodes |18]. In 
the five electrode planar trap design shown in Fig. [lib), the 
center and outermost electrodes are held at RF ground while 



FIG. 1 : Schematic illustrations of (a) a two-level linear RF Paul trap 
representative of what is currently used in quantum computation ex- 
periments, and (b) the five electrode planar ion trap suggested by 
Chiaverini et al. tv/h . Ions are trapped along the trap axis, shown 
as a dotted line. An RF potential is applied to the red electrodes to 
provide radial confinement, and DC potentials are applied to the blue 
control electrodes to provide axial confinement and to shuttle ions 
along the trap axis. Typical dimensions for current two-level traps 
are a slot width s of 200-400 jUm. 



the remaining two electrodes are biased with an RF potential 
for radial confinement. Either the center electrode or the out- 
ermost two electrodes can be segmented and DC biased for 
axial confinement. Such planar ion traps can be built using 
silicon VLSI technology and thus have the capability to scale 
to arbitrarily large and complex trap arrays 1 19]. 

Several challenges will have to be addressed, however, be- 
fore planar ion traps can be used for quantum computing. Per- 
haps the most significant challenge is that planar ion traps 
have trap depths that are only of order 1% that of a multi- 
level trap of comparable dimensions 1 17]. While this is not 
a problem once the ions are loaded and laser cooled to near 
the ground state of the motion, it makes loading from a ther- 
mal ion source at room temperature or above difficult. A sec- 
ond challenge is to show that it is possible to perform the 
three basic ion movement operations required to implement 
the Kielpinski et al. |3] ion trap quantum computer architec- 
ture in a planar trap: ballistic ion transport in straight lines, 
splitting and joining pairs of ions in a single trap, and shut- 
tling ions around corners in intersections. Further challenges 
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are that optical access for laser cooling is blocked along the 
axis orthogonal to the trap substrate, and that ion heating is 
enhanced first by the close proximity of the electrodes that is 
required to achieve a reasonable trap depth, and second by the 
more resistive semiconductor materials used for VLSI 1201 . 

In this work we address the first and second challenges. 
The latter challenges are addressed elsewhere Ellli. We 
investigate planar ion traps using printed circuit board (PCB) 
electrodes to trap macroscopic charged particles in a rough 
vacuum environment. We do not cool our ions, and imag- 
ing is accomplished by classical (off resonant) scattering of 
laser light. Because these experiments do not require an ultra 
high vacuum environment, we are able to achieve a fast cycle 
time for testing different electrode layouts. We optimize the 
electrode geometry to maximize the trap depth, and further 
demonstrate loading from a conventional four rod linear Paul 
trap acting as a reservoir of ions. Atomic ions could be laser 
cooled in the four rod trap before being transferred to the pla- 
nar trap to circumvent the problem of low trap depth. We test 
trap layouts with segmented center electrodes and segmented 
outer electrodes, and have accomplished all three of the move- 
ment operations required for the Kielpinski^t al. fs] architec- 
ture. While transport, splitting and joining |8], and movement 
through intersections 1 9] have already been achieved in multi- 
level traps, this work represents the first realization of these 
operations in planar ion traps. 

This paper proceeds as follows. Section goffers a general 
discussion of planar ion trap design and a detailed descrip- 
tion of our experimental setup. In Sec. |3l we characterize 
the macroscopic charged particles we use to test the traps. 
Section |4| presents the results of our investigation of planar 
ion traps including secular motion and ion movement exper- 
iments. In Sec.|5lwe discuss some of the various alternative 
trap geometries we have tested and propose a novel planar ion 
trap geometry that forms a lattice of point RF Paul traps. Fi- 
nally, in Sec. we summarize and suggest how to proceed in 
experimentally demonstrating atomic planar ion traps. 



2. EXPERIMENTAL DESIGN 

The design of a planar trap for trapping macroscopic par- 
ticles involves three main considerations: the ion source, the 
ion loading scheme, and the PCB design. Particle ion sources 
in common use include electrospray ionization 1I21I1 . laser ab- 
lation 1 22, 23], and piezoelectric particle generation |24]. The 
careful choice of an ion loading mechanism is important be- 
cause of the low trap depth of planar traps. The planar trap 
described in this paper is loaded using a four rod trap with a 
higher trap depth. Designing the PCB involves choosing the 
electrode dimensions to get the desired trap depth, ion height 
above the substrate, and secular frequencies; choosing the lay- 
out of the control electrodes to allow for performing the de- 
sired movement operations; and choosing a PCB substrate that 
has a high voltage breakdown threshold, does not accumulate 
charge easily, and is vacuum compatible in the case of a trap in 
a UHV environment. This section begins by presenting some 
general considerations for planar trap design and proceeds to 



describe the experimental setup. 

2.1. Trap design 

The dynamics of an ion in a linear RF Paul trap are deter- 
mined by solving the classical equations of motion ifisll . For 
an RF quadrupole electric potential 

2 2 

^{x,y,t) = ^^—^{U-Vcos{^t)) , (1) 
zrQ 

the X and y equations of motion take the form of Mathieu equa- 
tions: 

^ -^{a-2q cos(2t) ) jc = (2) 

and 

-4 - - 2qcos{2T))y = . (3) 
dt"^ 

Here t = Q,t/ 2 is a dimensionless time, the Mathieu param- 
eters q = {2QV)/{mrlQ?) and a = {4QU) / {mrlO?) are di- 
mensionless RF and DC voltages, and Q and m are the ion 
charge and mass. The motion is stable (i.e., the components 
of the ion position vector x and y are bounded in time) in re- 
gions of a — q parameter space. In particular, for <2 = the 
motion is stable for < ^ < qmax = 0.908. 

In the pseudopotential approximation, where ^ <C 1 , the ion 
motion along axis / can be decomposed into slow, large ampli- 
tude secular motion at the secular frequency cOt and fast, small 
amplitude micromotion at the RF drive frequency CI 12511 . For 
an arbitrary electric potential of the form 

= (t)RF{x,y,z)cos{at) ^(t)Dc{x,y,z) , (4) 

the secular motion is determined by a secular potential 

Ysec{x,y,z) = -^^^\V(l)RF{x,y,z)\^ ^Q(l)Dc{x,y,z) . (5) 

Using the secular potential we can calculate the height ro of 
the ion above the plane of the electrodes, the secular frequen- 
cies cOi in the harmonic region of the potential, and the trap 
depth y/secO- The trap depth is defined in the usual way as the 
energy an ion would need to escape the trap. The secular po- 
tential for the planar trap used in these experiments is shown 
in Fig. 121 

In order to determine the optimum relative sizes of the trap 
electrodes, we have calculated the secular potential numer- 
ically for the planar trap geometry defined in Fig. |2| over a 
range of dimensions using a two dimensional finite element 
electrostatic package named BELA |26]. Here Wc is the width 
of the center electrode, Wr is the width of the RF electrodes. 
Wo is the width of the outer electrodes, g is the width of the gap 
between the center and the RF electrodes, g' is the width of the 
gap between the RF and outer electrodes, h is the height of a 
planar electrode which is parallel to and above the trap elec- 
trodes, ro is the height of the ion above the trap electrodes. 
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FIG. 3: Trap depth versus electrode widths. The normaUzed trap 
depth do is plotted in color as a function of logio(wr/^o) and Wc/ro 
for h/ri = 706. The dashed lines are contours of constant g/rQ. Note 
that the maximum trap depth is obtained at Wc/ro ^0.5 and Wr/ro 
as large as possible. 



FIG. 2: Secular potential for a planar ion trap. The particular geom- 
etry depicted is that of the cross trap described in Sec. 12.21 The sec- 
ular potential is plotted on a linear color scale as a function of x and 
y, with blue representing the lowest and red representing the highest 
secular potential values. Secular potential values above twice the trap 
depth are truncated for clarity. The ion is located at the origin of the 
coordinate system. An RF potential Vcos{Q.t) is applied to the red 
electrodes, a DC potential U is applied to the blue top plate electrode 
(U = in this figure), and the black electrodes are grounded. 



and ri = {wc + w^) / 2 + g is a measure of the trap size. We set 
= g and Wo and varied w^, w^, and h. Note that h^oo 
for a true planar trap. Figures |3l and |4| show the normalized 
trap depth do defined by 



i2t/2 



Q'V 



do , 



and the normalized secular frequency defined by 

QV 



COi 



fi 



(6) 



(7) 



as functions of Wc/ro and Wy/ro for h/r\ = 706. Note that 
do= fi = ^ for a perfect quadrupole trap and that we consider 
the secular frequency in the harmonic region of the potential 
where it is the same along both axes (this is appropriate when 
the ions are at low temperature). 

Using Figs. [3land|4l we can design a planar trap to max- 
imize the trap depth subject to the experimental constraints. 
As an example, consider a planar ion trap for ^^Sr+. The ion 
height ro is set by either the acceptable ion heating rate or 
by how tightly the cooling and other lasers can be focused 
such that they do not scatter off the surface of the trap. Sup- 
pose we pick ro = 500 /im. The electrode widths Wc and Wr 
should be chosen to maximize the trap depth using Fig. [3j 
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FIG. 4: Secular frequency versus electrode widths. The secular fre- 
quency at zero ion temperature fi is plotted in color as a function of 
\ogio{wr/ro) and Wc/ro for h/ri = 706. The dashed lines are con- 
tours of constant g/ro. 



subject to the constraint that the gap between them g must be 
large enough that the RF voltage does not induce electrical 
breakdown. In general, the optimum center electrode width is 
around 0.5ro, and the optimum RF electrode width is as wide 
as possible given practical considerations such as the capaci- 
tance of a large electrode. For our example, suppose further 
that g O.Sro 250 jim and Wr = ro = 500 /im. Figs. [H 
andil then, give Wc = 0.58ro = 290 fim, do = 0.0065, and 
fi = 0.28. At y = 500 V and ^1 = 2;r X 10 MHz, this corre- 
sponds to Xj/secO = 0-47 eV and cO/ = 27r x 1.1 MHz. 

We have found that by adding a planar top electrode (i.e., 
choosing h to be finite) biased at a positive DC potential U we 
can increase the trap depth by up to a factor of a 40 for a given 
RF amplitude and frequency. Figure |5] shows the normalized 
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trap depth defined by 



Vsec.O ^ 



g2y2 
o — 

Amr\Q? 



(8) 



as a function of the normaUzed voltage on the top plate, 

-U , (9) 



u = 



Amr\Q?' r\ 



QV^ h 



for Wc/ri =Wr/ri =0.6 and several values of /z/ri . Note that 
the trap depth is normalized differently for this plot than it 
was for Fig.[3l This is because ro is a logical starting point for 
selecting the electrode dimensions, but r\ is the fixed length 
scale once the trap is built. At w = 0, the smallest barrier in 
the secular potential over which the ion can escape is in the 
positive y direction. As u increases, the minimum energy es- 
cape path shifts first to the sides and then to the negative y 
direction. For a given geometry, the trap depth increases with 
u until the minimum energy escape path is straight down to 
the center electrode, then decreases with a further increase in 
u. The maximum trap depth occurs at w ^ 1 and is given by 
d\ 0.29, or a little more than one fourth the trap depth of a 
standard four rod ion trap with similar dimensions. The secu- 
lar potential depends only weakly onh/r\ because the effect 
of nonzero u is well approximated by a constant electric field 
in the regime h/r\ ^ 1 that we consider. Note that because 
the DC potential pushes the ion off of the RF quadrupole null, 
there will be increased micromotion. To first order in a har- 
monic potential, the micromotion amplitude is 



(10) 



where is the displacement of the ion from the RF null. In 
practice we envision turning the DC potential on for loading, 
then gradually turning it off as the ions are laser cooled. Con- 
tinuing our previous example, with a top plate electrode h=\ 
cm above the trap biased at = 100 V, we get a trap depth 
Wsec.Q = 4.6 eV and a micromotion amplitude of 15 /im using 
Eq. dTol . A direct numerical calculation of the classical ion 
trajectory along the j-axis predicts a micromotion amplitude 
of 20 /im, confirming that Eq. dTol is a reasonable approxima- 
tion despite the rather non-harmonic nature of the potential. 



2.2. Experimental Apparatus 

In this section, we discuss the specifics of the trap's design 
and construction, our source of ions, the method for loading 
the planar trap, the environment used for performing our ex- 
periments, and our methods for measurement and control of 
ions. 

The planar trap we investigated is made up of four straight 
arms joined at a cross intersection (Fig.|6|). It is a printed cir- 
cuit board (PCB), made using standard techniques. The elec- 
trodes are tin coated copper, and the substrate is GML-1000, a 
microwave laminate. The trap is a five electrode design, with 
outer electrodes segmented to control the axial potential. Con- 
nections to these electrodes are made on the bottom layer of 
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FIG. 5: Trap depth versus top plate location and potential. The 
normalized trap depth di is plotted as a function of u for Wc/n = 
Wr/ri = 0.6 and several values of h/ri. For u < 0.2, the ion escapes 
straight up in the positive y direction; for 0.2 < w < 1, the ion es- 
capes out to the side; and for 1 < w the ion escapes straight down in 
the negative y direction. Note that the maximum trap depth di ^ 0.29 
is obtained at w ?^ 1 . 



the PCB using surface mount headers. Most opposing pairs 
of control electrodes are electrically connected, but the two 
pairs nearest the intersection on each arm are electrically in- 
dependent for finer control. The center electrodes are also seg- 
mented at the intersection to provide additional control in that 
region. The middle three electrodes are all 1.27 mm wide, the 
electrode spacing is 0.89 mm, and the outer control electrodes 
are 2.5 mm long. We do not use the optimal electrode sizes 
as discussed in Sec. l2.1l due to the practical considerations of 
printed circuit board manufacture, specifically minimum fea- 
ture and drill sizes. 

Electrospray ionization (ESI) provides ions to the trap. This 
technique, commonly used with linear quadrupole filters for 
mass spectroscopy, applies high voltage to a liquid solution 
at a sharp tip. Strong electric fields at the tip blow off fine 
droplets of solution, and as solvent evaporates from these 
charged droplets, self-repulsion breaks them up into smaller 
particles, eventually producing individual charged solute par- 
ticles. Following 1271] . we begin with a 5% suspension of 0.44 
/im diameter aminopolystyrene spheres. We then prepare a 
solution buffered to pH = 3.9, add methanol to produce a 4:1 
methanol/solution mixture, and add the spheres to produce a 
0.05% suspension. The suspension is then placed in a sealed 
bottle, and pressurized to 50-100 kPa. This forces the liquid 
through a 0.45 fim filter to block any clusters of microspheres, 
and out the electrospray tip. The tip itself is made from fused 
silica capillary tubing, which is heated, stretched to produce 
a neck, and cleaved to produce a 100 jim opening. We apply 
4 kV directly to the liquid with a copper wire inserted in the 
fluid near the tip (Fig.Q. The spray plume from the tip is di- 
rected through a grounded mask and into the end of a four rod 
trap. 

By arranging the rods so that they match the strips of a pla- 
nar trap, we interface the four rod trap directly to the planar 
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FIG. 6: Top view of the trap. Opposing electrodes in straight sections 
of the trap are electrically connected, but near the intersection both 
the outer and center electrodes are separated to provide finer control. 
Electrical connections are made via surface mount headers on the 
underside. The RF loop at each end helps prevent ions from leaking 
out axially. 
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FIG. 7: The ESI apparatus. Microspheres suspended in a mixture of 
methanol and water are forced from a pressurized bottle, through a 
0.45 /im filter and past a copper wire with +4 kV applied voltage, 
and out the capillary tip. The charged spray passes through a mask 
and into a four rod trap. 




FIG. 8: Four rod to planar interface. A slot cut in the planar trap 
allows for a common trap axis. The large particles form a linear 
Wigner crystal immediately and are pushed into the planar trap. 




Laser Axis 2 



FIG. 9: Top view of the vacuum chamber. Ions are loaded into the 
left side, and then a flange is screwed into place and the air evacuated 
slowly. Two green lasers illuminate the trap axes. 



trap. A slot cut in the side of the PCB allows us to lower the 
four rod trap to a point where its trap axis coincides with that 
of the planar trap (Fig.[8l). Mutual repulsion of ions loaded in 
the far end of the four-rod trap forces them along, eventually 
pushing them into the planar trap. When completely full, the 
planar trap can hold ~50 ions. 

Electrospray ionization must be performed at or near atmo- 
spheric pressure, and experiments on free trapped particles re- 
quire a vacuum, so we enclose both four rod and planar traps 
in a custom built clear acrylic box (Fig.|9ll, and load through 
an open side of the box. We perform all loading at atmo- 
spheric pressure, then seal the box shut and carefully pump 
down. We use a needle valve to control flow rate during pump- 
ing, since any significant air flow can push ions out of the trap. 
The lowest attainable pressure in this enclosure is about 15 
Pa (10~^ torr). Green (532 nm) laser light directed along the 
trap axes illuminates ions by classical (off resonant) scatter- 



ing. The ions are clearly visible by eye, and we photograph 
them using an ordinary camera lens mounted on a CCD cam- 
era. 

Also shown in Fig.|9lis a metal plate installed above the trap 
to increase its trap depth and shield from stray static charge. 
Machined slots in this plate allow us to view ions from above, 
but we can also see them from the side. This slotted top plate 
design has the advantage of masking laser scatter, but for more 
complex trap topologies it might be easier to use a transparent 
conductor such as indium tin oxide or a thin film of gold |28, 
l29[ I3QIi . This could be deposited on a glass plate or directly 
on the vacuum window or an imaging optic inside the vacuum 
chamber. 

The electronics used to generate the RF and DC potentials 
for the trap are shown in Fig. [lO| Frequencies for which the 
microspheres are stable range from a few hundred hertz to sev- 
eral kilohertz, depending on the pressure. This is in the audio 
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FIG. 10: Control electronics for the planar trap. 



frequency (AF) range, but we will continue to refer to it as RF. 
Since these frequencies are easily synthesized, we use a func- 
tion generator as the source, then use an active high voltage 
amplifier to reach the target voltage range. A small tickle sig- 
nal can be added to the main RF signal before amplification 
to probe for resonances in ion motion, allowing us to deter- 
mine secular frequencies (Sec. 14.21) . In this design, the upper 
limit to RF voltage is set by arcing between RF and ground 
electrodes, which occurs at about 400 V amplitude. Typical 
operating conditions are 250 V and 1.5 kHz. 

To provide DC control of the trap's 48 independent elec- 
trodes, we use a board with 24 digital output channels to ad- 
dress six serial octal digital to analog converter (DAC) inte- 
grated circuits. The outputs of these DAC chips are then am- 
plified to ±20 V, sufficient to axially confine or to move these 
particles. 



3. ION CHARACTERIZATION 



18 
16 
14 

)12- 
!lO 



CD 

E 



I 



Ua 



Charge/mass (e/u) 



X 10" 



FIG. 1 1 : Charge-to-mas s distribution of ESI-generated ions. The 
mean 2/misl.l7x 10~^ e/sLum and the standard deviation is 0.49 x 
10"^ e/mnu. 



= 2n x2 kHz and V = 250 V. With the chamber pressure 
at 70 Pa, the RF frequency was lowered while observing the 
ejection frequency Q^ej of each ion. The charge-to-mass ra- 
tio was calculated using the formula Q/m = qmaxi"o^lj/ 
where q^ax = 0.908, the value of the Mathieu q parameter 
at the boundary of the first stability region of a quadrupole 
trap. Figure [n] shows a spectrum obtained for 178 ions with 
mean 2/m = 1.17 x 10~^^/amu (electronic charges per nu- 
cleon mass). 



Our goal in this work is to study the properties of planar 
traps rather than the ions loaded into them. To be able to draw 
conclusions about atomic ion traps from work on macroscopic 
ions, we must (i) scale away the charge-to-mass dependence 
from measured parameters and (ii) understand the effect of 
damping on macroscopic ions since it is negligible for atomic 
ions. Section lTTl discusses measurement of the 2/m spectrum 
of the ions. Section lT2l discusses the effects of air drag on the 
ions' stability and motion. 



3.1. Charge to mass ratio 

The principal drawback of using macroscopic charged par- 
ticles to investigate a trap design is that, unlike atomic ions, 
these particles are not identical. Electrospray ionization (ESI) 
is a soft ionization technique that produces multiply charged 
ions with a significant spread in the charge state. The varia- 
tion in the mass of the microspheres 10%) also contributes 
to broadening the charge-to-mass distribution. Nevertheless, 
we are able to keep the 2/m spectrum constant by maintain- 
ing constant electrospray pressure, voltage, capillary diameter 
and loading parameters. 

To measure the 2/m spectrum of the ions, groups of about 
five ions at a time were loaded into the four rod trap with 



3.2. Damping of ion motion by background gas 

When trapping macroscopic ions at atmospheric pressure 
or in rough vacuum, the background gas exerts a drag force 
on the ions that both stabilizes the radial motion and slows 
down the axial ion movement operations. It is important to 
verify that the background pressure in our experiments is low 
enough to make accurate assessments of the trap performance. 

For small Reynolds number Nrc = 2pxR/il and small 
Knudsen number Nku = A/7?, the drag on a sphere is given 
by Stokes' law 

FD = -67riiRx. (11) 

Here p is the density of the gas, R is the radius of the sphere, 
II is the dynamic viscosity, A is the mean free path, and x is 
the ion velocity along the x axis. When the Knudsen number 
is large. Stokes' law overestimates the drag because the flow 
is no longer continuous. We divide Eq. JTTl by an empirical 
slip correction factor 

C{NKn) = l^NKn (l . 165 + 0.483^-0-99^/^^") (12) 

to obtain an expression for the drag force which is valid at 
large Knudsen number L3L] . 
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FIG. 12: Drag coefficients describing ion stability and the speed of 
linear ion shuttling versus background pressure. The dimensionless 
drag coefficients b and c for stability and ion shuttling, respectively, 
are plotted as functions of the vacuum chamber pressure on a log- 
log scale. The drag coefficients are computed using R = 0.22 /im, 
2 = 5.3 X 10-^^ C, m = 4.7 x 10"^^ kg, ji = 1.83 x 10"^ kg m-^ 
s-\ A = (67.3 nm)(105 Pa//?), a = In x 5 kHz, d = \ mm, and 
E = (1 V)/(l mm); which is appropriate for the cross trap assuming 
that we are trapping single microspheres. The circles indicate the 
values of the drag coefficients at the pressure used in our experiment. 

We determine the effect of drag on stability by considering 
a modified dimensionless equation of motion 

d^x dx 

^b—- + (a - 2^cos(2t))x = (13) 
where the dimensionless drag coefficient 



C{NKn)m^' 

A similar equation holds for the y motion. For /? = 0, this 
reduces to Eq. (|2|). Numerical computations verify that the 
region of stable trapping ma — q parameter space grows with 
increasing b |32,l33l- In particular, the maximum stable Math- 
ieu q parameter at <2 = goes from qmax = 0.908 at /? = to 
^max = 1.05 at /? = 0.45. The drag parameter b is plotted as a 
function of vacuum pressure p in Fig.[l2| At 70 Pa (0.5 torr), 
b = 0.45 is already less than 1 so air drag has only a small 
effect on the trap stability in these experiments. 

The speed of ion motion in our shuttling experiments, how- 
ever, is significantly affected by drag. The physical reason for 
this is that the characteristic time for shuttling experiments is 
much longer than the characteristic time for stability (1/^^). 
As a simple model, consider an ion which starts from rest at 
z = in an axial potential 

, . . / E(d-z) , z<d 
^W = \ , z>J ^^^^ 

and suppose we measure the ion velocity when z = d. In this 
model the length scale d is of the order of the axial dimension 



of a control electrode, and the electric field E is of the order 
of the electric potential applied to a control electrode divided 
by its axial dimension. Drag is important when the dimen- 
sionless drag parameter c = y-fj > 1 where 7= — F/)/(mz) and 
is the time to reach z = d. Here is the solution to the 
transcendental equation 

QE ( e-^d - 1 \ 

— Uj + . (16) 

ym \ 7 J 

The dimensionless drag coefficient c is plotted as a function of 
pressure in Fig.[l2l At 70 Pa, c = 4.9 so drag plays an impor- 
tant role in ion shuttling experiments in the cross trap despite 
the fact that the trap stability is not affected by damping. 

For large damping (c 1), z{z = d) approaches the 
terminal velocity {QE)/{ym). The drag coefficient 
\/C{NKn) ~ I/Nku ~ SO z{z = d) ^ I /p. This is experi- 
mentally verified in Sec. 14.31 



4. PLANAR TRAP PERFORMANCE 

When damping forces can be neglected, the Mathieu equa- 
tions of motion of a trapped ion depend only on dimension- 
less parameters. It follows that the dynamics of macroscopic 
charged particles, viewed in the appropriate time scale (a mi- 
cromotion period), are identical to that of atomic ions, and 
can be fully explored without the much more demanding ex- 
perimental requirements of trapping atomic ions |34]. We 
have found macroscopic charged particle planar traps to be 
a rapid and accurate test bed for investigating traps of differ- 
ent geometries (e.g. three electrode, five electrode) and con- 
trol electrode layouts (e.g. segmented center, segmented outer 
electrodes). In this section, we present the results of an ex- 
perimental investigation of planar trap performance using this 
test bed. 



4.1. Ion height above trap substrate 

A planar trap with a segmented center electrode can be used 
to control the height of individual ions above the substrate by 
applying DC voltages to the electrode segments. One advan- 
tage of such control in a quantum computer is the ability to 
perform single qubit gates on different ions using a single sta- 
tionary laser beam. The laser beam would have to be parallel 
to the ion chain but slightly raised or lowered. By changing 
the height of individual ions, they can be brought into the path 
of the beam to perform single qubit gates. 

We measured the height of an ion above the trap sub- 
strate using a CCD camera mounted on a calibrated translation 
stage. The micromotion of the ion in the y direction causes 
the images (which have an exposure time much longer than 
the period of the RF drive) to show a streak where the ion is 
located, so we actually measured the positions of the top and 
bottom of the ion motion. Figure [21 shows the average of the 
positions of the top and bottom of the ion motion as a func- 
tion of the DC potential on the center electrode. We measured 
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FIG. 13: Ion height above the trap substrate, ro, is plotted as a func- 
tion of the DC potential on the center electrode. The solid line is 
the expected ion height from numerical computations of the secular 
potential as described in Sec. l2.1l using the measured charge-to-mas s 
ratio of the ion. 



4.2. Secular frequencies and trap geometric factors 

We determine the secular frequency of an ion by adding a 
Fourier component of frequency Qt and amplitude Vt to the 
RF signal. Typically, Vt is between 0.5 V and 2 V. By sweep- 
ing Q,t^ we are able to observe a sharp increase in the size 
of the ion trajectory along the direction of the secular mode 
excited by the drive signal as seen in Fig. [21 Although the 
drive signal is applied to the RF rods, indirect coupling of 
the drive to the axial motion of the ion is sufficient to excite 
axial modes. The exact axial secular frequency depends on 
which trap electrodes are used as endcaps and what voltage 
is applied to them, but it is typically between 15 Hz and 30 
Hz. The transverse secular frequencies are of more interest 
because they depend directly on the trap geometry: the geo- 
metric factor fi/rl in Eq. depends only on electrode di- 
mensions and their spacing. We were not able to resolve sep- 
arate resonant frequencies for cOj and cOy, which indicates that 
the ion energy is much less than the trap depth. Figure fT6l 
shows the dependence of the x secular frequency on V and Q,. 
Both dependences behave as expected from the pseudopoten- 
tial approximation. 





FIG. 15: (a) Top view of an ion trajectory from secular resonances 
with V = 250 V and Q. = 2n x2 kHz. Only the center ground elec- 
trode is seen in these pictures because the grounded top plate elec- 
trode masks the other electrodes from view. The width of the ground 
electrode is 1.27 mm. (b) Transverse mode excited at a drive fre- 
quency of 288 Hz. (c) Axial mode excited at 25 Hz. 
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FIG. 14: Ion motion amplitude. The peak to peak amplitude of the 
ion motion in the y direction is plotted as a function of the DC poten- 
tial on the center electrode. 



350 = 
^00 



g/m for the ion by lowering the RF drive frequency until the 
ion became unstable as described in Sec. 13.11 The expected 
ion height shown in Fig. [21 is calculated using this value of 
Q/m and a numerical computation of the secular potential. 

Figure[l4| shows the amplitude of the ion motion as a func- 
tion of the DC potential on the center electrode. The micro- 
motion is minimized when the ion is located at the null of the 
RF quadrupole electric field, which occurs when the potential 
on the center electrode is zero. Even without a bias on the cen- 
ter electrode, however, there is still a substantial amplitude of 
ion motion because we do not cool the secular motion of our 
ions. 
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FIG. 16: Measurements of the secular frequency versus (a) RF am- 
plitude at a fixed RF frequency of 1.5 kHz (on a linear scale) and (b) 
RF frequency at a fixed RF amplitude of 250 V (on a log-log scale). 
The solid lines are fits to Eq. with ^fi as a fit parameter. 

To confirm the accuracy of our simulations of pseudopo- 
tentials in planar traps, we measured the ratio of the geomet- 
ric factor fx/rl in the planar trap to that of the loading four 
rod trap. The large spread in the Q/m distribution does not 
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FIG. 17: (a) A strobe image showing shuttling of a single ion, taken 
by modulating the current to the illuminating laser. The laser is pe- 
riodically turned on for 2 ms and off for 5 ms. The ion is moving 
from left to right; there is a short acceleration period after which the 
ion reaches a maximal velocity and then decelerates. This shuttle 
was performed by applying a potential of 5 V to the nearest control 
electrode that is to the left of the ion. (b) Distance moved by the ion 
versus elapsed time, fitted with a spline to guide the eye. 

allow for the determination of the geometric factor for each 
trap separately. However, we can measure the ratio of the ge- 
ometric factors of the traps by first measuring (0^ of an ion 
in the planar trap, then applying a sequence of voltages on 
the control electrodes to move that ion into the four rod trap, 
where COx is measured again. By doing several experiments 
of this kind, we determined the geometric factor ratio to be 
0.40 ±0.01. Numerical simulations of this setup were also 
performed using BELA, giving for this ratio the value 0.404, 
which is within the error bar of the experimental result. 

4.3. Ion movement 

Three important ion movement primitives for a multiplexed 
ion trap quantum computer are shuttling, turning corners 
and splitting/joining two ions |3]. We have demonstrated 
these three operations in planar traps, using our macroscopic 
charged particle system. 

Figure [n] shows a strobe image of shuttling a microsphere 
in the planar trap at 40 Pa (0.34 torr). The shuttling shown 
was performed by raising the potential of the nearest control 
electrode to the ion. In this case, the speed of the ion and its fi- 
nal position are determined both by the electrode potential and 
the drag force on the ion. Figure[^ shows how the maximum 
speed the ion achieves during its motion depends on these two 
parameters. 

We have also shuttled ions using another scheme where the 
nearest control electrodes on both sides of an ion are used to 
create a confining well along the axial direction. By moving 
this well along the trap, we had more control over the acceler- 
ation, speed, and final position of the ion. 

Splitting a pair of ions was performed by using the con- 
trol electrodes to introduce a potential hill that pushes the ions 
apart. The ion spacing is approximately equal to the axial 
dimension of one control electrode, so this technique is sat- 
isfactory. Typically, an electrode between the ions is raised 
to 5 V for separation and lowered back to V to join them 
again. For performing two qubit gates with atomic ions, it is 



100 




5 10 15 20 

Electrode voltage (V) 



10° 

Chamber pressure (torr) 



FIG. 18: The maximum speed Vmax an ion attains during a shut- 
tling operation versus (a) the potential Uc applied to the nearest con- 
trol electrode at constant pressure (40 Pa) fit to a 3/4 power law 
(vmax =Aul^^ where A = 9.55 cm s~^V~5) (b) the chamber pres- 
sure at constant control electrode voltage (5 V). Ion speed is found to 
be inversely proportional to the chamber pressure. This behavior is 
expected when (1) the mean free path is much greater than the radius 
of the particles and (2) ytac ^ 1 where y the damping factor and tac 
is the acceleration time of the ion. When the damping is low enough 
such that (2) no longer holds, the ion speed is expected to reach a 
pressure-independent value. 



desirable to bring the ions as close to each other as possible, 
because gates are fastest when the frequencies of the two ax- 
ial normal modes of the ions are well separated |35]. In that 
case, the ion spacing is likely to be smaller than the electrode 
spacing, and it has been shown that it is advantageous to use 
the DC electrodes to produce an electric octopole moment for 
fast separation |36]. 

Figure [191 shows ion turning a corner in the planar trap. 
Simulations of the pseudopotential indicate that there are po- 
tential hills in the RF nodal lines in each of the four arms close 
to the intersection, with a point node right at the intersection. 
To turn a corner, it is necessary for an ion to overcome the 
hills in both the source and destination arms. We have used 
the following pulse sequence to do that: (1) lower the center 
electrode of the destination arm to -2 V, then (2) raise a con- 
trol electrode in the source arm to 5 V to push the ion over 
the hill in the source arm. We find that the ion has enough 
kinetic energy after traversing the first potential barrier to also 
overcome the second one. 
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(a) (b) 

FIG. 19: Turning a corner. The arrows point to the position of the 
ion (a) before and (b) after turning the corner. This operation takes 
about 50 ms. 



5. VARIATIONS ON TRAP DESIGN 

Macroscopic charged particles and PCB traps allow us to 
rapidly test a variety of planar trap layouts. We have focused 
on the loading and transport properties of multi-zone linear 
traps that are joined at intersections. In addition to the cross 
trap described in this paper, designs under investigation have 
as many as 100 zones and include corners, three way junc- 
tions, and four way junctions. These elements are necessary 
to generate the ion trap geometries that arise in complex quan- 
tum computation geometries |3], particularly those involved 
in realizing scalable, fault-tolerant quantum computation cir- 
cuits (371]. 

In addition to fault-tolerant quantum computation, ion traps 
are promising candidates for quantum simulations. Porras 
and Cirac have proposed using the motional modes of cou- 
pled ions to simulate Bose-Hubbard models |38]. The planar 
version of a point Paul trap is a natural way to implement this 
scheme for a two dimensional system. Starting with a plane of 
RF voltage one can place DC electrodes at arbitrary positions 
defining a planar point trap at each position. The distance 
between trapping points controls the strength of vibrational 
coupling and the layout of the electrodes can be used to create 
both ordered and disordered systems. 

Using this approach, we have implemented a hexagonal lat- 
tice of traps shown in Fig.|20l this lattice successfully trapped 
macroscopic ions. For experiments with single atom ions, 
the lattice spacing needs to be ~50 /im to be near interest- 
ing phase transitions in condensed matter systems such as 
that studied by Porras and Cirac. Such a trap can be oper- 
ated as described above with the ions trapped at the lattice 
points as shown in Fig. |20| If the RF and ground electrodes 
are switched, the ions are trapped along the honeycomb lines, 
and applying an offset voltage to the RF on the dots moves the 
ions through the lattice. This increased coordination between 
ion positions may improve fault- tolerant thresholds. 



6. CONCLUSIONS 

We have experimentally demonstrated the planar ion trap 
design using macroscopic charged particles as a rapid devel- 




FIG. 20: Top view of a two dimensional Paul trap array. Ions can 
either be confined to points or free to move along lines, depending on 
the RF and ground connections. The connections shown (top inset) 
trap ions above the dots. Bottom inset: Microspheres are trapped in 
the 2D array, and seven are shown illuminated by the laser beam. In 
contrast with the cross trap, the four rod loading stage is out-of-plane, 
and is visible above and to the right. 

opment test bed. We have addressed the challenges of ion 
loading and ion motion given the smaller relative trap depth 
of planar traps. 

The challenge of loading ions into the planar trap was 
solved by coupling a traditional four rod trap to a planar trap. 
Ions at high kinetic energy were first loaded into the linear 
trap and then offloaded onto the planar trap. This method al- 
lows us to load the planar trap with the efficiency of loading a 
traditional linear trap. 

Furthermore, we have found that the addition of a charged 
conductive plane above the trap can be used to increase the 
trap depth. In the atomic case the ions will be initally loaded 
at a position of high micromotion limiting the effect of laser 
cooling. In our design example this micromotion has an am- 
plitude of 0.04ro when the potential on the top plane is such 
that the trap depth is increased by a factor of 10. However, 
after initial trapping the micromotion can be minimized and 
the cooling maximized simply by grounding the conductive 
plane. 

Controlled movement of the ions in the planar trap was 
accomplished despite the small trap depth. Using increas- 
ingly complex traps, we have been able to perform all the 
fundamental movement operations required in a multi-zone 
architecture: splitting and joining ion chains and moving ions 
around corners and through four way intersections. Before 
this work, ion movement through a four way intersection was 
predicted to not be possible. 

Questions remain about the control of ions at low vacuum. 
Our experiments were performed at a vacuum where the back- 
ground pressure no longer effects trap stability. However, the 
pressure still contributed a significant drag term to the linear 
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motion. 

Additionally, we have demonstrated that planar traps can 
be used to produce a wide range of geometries. As an exam- 
ple, we have trapped in a hexagonal lattice. A similar trap for 
atomic ions could be used to simulate two dimensional quan- 
tum simulations. One can then controllably introduce disorder 
into the system by application of voltages or changes in trap 
fabrication. 



Planar ion traps are a general tool with many applicatuions 
in quantum simulation, quantum computation, and mass spec- 
trometry. Using modern two dimensional fabrication tech- 
niques, the traps described here can be reconstructed to be 
compatible with UHV and atomic ions. Open questions about 
the effect of the surface on ion heating can then be addressed. 
Planar traps or "ion chips" offer an exciting route to study new 
atomic ion physics. 
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